Atmos. Chem. Phys. Discuss., 9, 27465-27483, 2009 
www.atmos-chem-phys-discuss.net/9/27465/2009/ 

© Author(s) 2009. This work is distributed under 
the Creative Commons Attribution 3.0 License. 



Atmospheric 
Chemistry 
and Physics 
Discussions 


This discussion paper is/has been under review for the journal Atmospheric Chemistry 
and Physics (ACP). Please refer to the corresponding final paper in ACP if available. 


Aerosol-cloud interaction determined by 
both in situ and satellite data over a 
northern high-latitude site 

H. Lihavainen 1 , V.-M. Kerminen 1 , and L. A. Remer 2 

1 Finnish Meteorological Institute, P.O. Box 503, 00101, Helsinki, Finland 
2 NASA Goddard Space Flight Center, Greenbelt, MD, 20771 , USA 

Received: 4 November 2009 - Accepted: 27 November 2009 - Published: 21 December 2009 
Correspondence to: H. Lihavainen (heikki.lihavainen@fmi.fi) 

Published by Copernicus Publications on behalf of the European Geosciences Union. 


27465 


Abstract 

The first aerosol indirect effect over a clean, northern high-latitude site was investigated 
by determining the aerosol cloud interaction (ACI) using three different approaches; 
ground-based in situ measurements, combined ground-based in situ measurements 
5 and satellite retrievals and using only satellite retrievals. The obtained values of ACI 
were highest for in situ ground-based data, clearly lower for combined ground-based 
and satellite data, and lowest for data relying solely on satellite retrievals. One of the 
key findings of this study was the high sensitivity of ACI to the definition of the aerosol 
burden. We showed that at least a part of the variability in ACI can be explained by 
io how different investigators have related different cloud properties to “aerosol burden”. 


1 Introduction 

Atmospheric aerosol particles influence clouds, and thereby the Earth’s climate, by 
altering the albedo, lifetime and precipitation patterns of clouds (e.g. Twomey, 1974; 
Lohmann and Feichter, 2005; Stevens and Feingold, 2009). Collectively, these effects 
is are called aerosol indirect effects. The first aerosol indirect refers to the situation in 
which a cloud liquid path is held constant and additional particles are introduced into 
the cloud, resulting in more numerous but smaller cloud droplets and a more reflective 
cloud. 

The magnitude of the first indirect effect in the global atmosphere is poorly quantified 
20 (Lohman and Feichter, 2005; Penner et al., 2006; Quaas et al., 2008). One of the 
main reasons for this is that the response of cloud droplet number concentration to 
aerosol perturbations depends in a complicated manner on both aerosol properties 
and prevailing meteorological conditions (e.g. Shao and Liu, 2009). Experimentally, the 
first indirect effect can be determined most easily by using the concept of aerosol cloud 
25 interaction (ACI) that relates the cloud optical depth, droplet effective radius or droplet 
number concentration to the aerosol burden (Feingold et al., 2001). The magnitude 
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of ACI can be determined based on either in situ or remote sensing data, or on any 
combination of these two. 

In this study, we will investigate the first indirect effect over a clean, northern high- 
latitude site. By determining the aerosol cloud interaction (ACI) using three different 
5 approaches that rely on ground-based in situ measurements and satellite retrievals, 
we are searching for an answer to the following questions: (1) can we detect the first 
indirect effect by all these methods at our site?, (2) do these methods produce similar 
values for ACI, and if not, why?, (3) is the value of ACI affected by the definition of 
aerosol burden, and (4) what implications our findings might have on future studies on 
io this subject. 


2 Materials and methods 

2.1 Measurements 

The data analyzed here include different data sets and approaches. Some of these 
have been published elsewhere, so only a brief description is given here. A general 
is overview of the ground site and typical conditions there has been given by Hatakka et 
al. (2003). 

2.1.1 Ground based in-situ measurements 

The data set 1, called “Cloud cases", is based on measurements conducted between 
2000 and 2002 at the Pallas area in northern Finland (Komppula et al., 2005). The 
20 approach takes advantage of two stations with horizontal distance of 6 km at different 
elevations. The higher-altitude station, Sammaltunturi, is inside a cloud about 5% of 
the time, whereas the lower-altitude station, Matorova, is always below the cloud layer. 
There are two similar Differential Mobility Particle Sizer (DMPS) measuring particles in 
a dry size range 7-500 nm with 30 discrete size bins. The size distribution analyzed 
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here are for individual cloud cases lasting from a few hours to a few days. The mea- 
surement systems and data analyses have been described in detail by Komppula et 
al. (2005). 

The data set 2, “First PaCE’ (The First Pallas Cloud Experiment), is based on an 
5 intensive measurement campaign carried out in the fall of 2004 in the Pallas area. The 
basic measurement approach is the same as in the data set “Cloud cases”, but there 
were some additional measurements at both stations. For example, the cloud droplet 
number size disteribution was measured directly with a Forward Scattering Spectrom- 
eter Probe (FSSP) in Sammaltunturi. Data from the DMPS systems analyzed here are 
io 1-h averages. More details on the measurement systems and data analyses can be 
found in Lihavainen et al. (2008). 

The data set 3, “Second PaCE” (The Second Pallas Cloud Experiment), is based on 
an intensive measurement campaign conducted in the fall of 2005 in the Pallas area. 
The measurement approach is different from that during the “First PaCE” described 
is above. The interstitial and total aerosol properties were both measured at the Sam- 
maltunturi station using different inlet systems. The total air inlet system evaporates 
cloud droplets by heating. Interstitial aerosol properties were measured from a sample 
line that has a PM 2 5 inlet nozzle. Similar DMPS systems were attached to both lines. 
Additional measurements included those by an aerosol mass spectrometer and FSSP. 
20 The DMPS data analyzed here are one-hour averages. More details on the measure- 
ment systems and data analyses can be found in Anttila et al. (2009) and Kivekas et 
al. (2009). 

The data set 4, “2007 Data", was measured during a few months beginning in Au- 
gust 2007. The measurement system was similar to that during the “Second PaCE”. 
25 The only additional measurement device was a FSSP that was running on the side with 
the DMPS systems. DMPS data analyzed here are one-hour averages. 
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2.1 .2 Remote sensing and ground base measurements 

For this data set (number 5), we took advantage of combining the long-term ground- 
based measurement at Sammaltunturi and cloud retrievals from the MODerate resolu- 
tion Imaging Spectroradiometer (MODIS) on both the Aqua and Terra platforms. The 
s data analyzed in this context covers the years 2000-2005. The cloud properties used 
here are an average of the 1 kmxl km retrieved effective radius and optical thickness 
(Platnick et al., 2003) of a 5 kmx5 km box centered on the Sammaltunturi station. The 
data from MODIS were combined with the closest one-hour average data measured at 
Sammaltunturi, provided that it met the following requirements: 

io (a) The cloud was a low-level or boundary-layer cloud, with the limiting cloud top 
pressure set to 780 hPa. 

(b) The liquid water path (LWP) should be constant when using the cloud droplet 
effective radius (r e ) or cloud optical thickness (COT) in defining the aerosol-cloud 
interaction (ACI, see Eqs 1 and 2 below). This was not possible due to the small 

is number of data, so we decided to use a limited range of LWP in order to get as 
much data as possible. The chosen arbitrary range was from 100 to 200 gm" 2 . 

(c) The Sammaltunturi station was below the cloud base. This information we got 
from the visibility sensor which measures the visibility in the range from 10 m to 
50 km. The limiting value for one-hour average visibility was set to 1 km. 

20 (d) Since the ground-based data sets 1-4 are mostly from late summer and fall, the 

analysis here was restricted to days of year (doi) larger than 140. 

(e) No rain was observed at Sammaltunturi. 

2.1 .3 Remote sensing 

For this data set (number 6 ), we expanded the analysis from the local stations and 
25 took again advantage of both MODIS aerosol (Remer et al., 2005; Levy et al., 2007) 
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and cloud (Platnick et al. 2003) retrievals from both Terra and Aqua platforms. The 
aerosol and cloud properties were averaged over a 1 °x 1 ° grid. The area for which 
the calculations were done was from 65° to 70° N and from 20° to 30° E. The time 
period was from 2003 to 2005. The same limiting values were applied for the cloud top 
5 pressure, seasons and LWP as for the data set 5. 

2.2 Determination of the aerosol-cloud interaction 


The first aerosol indirect effect is defined as the change in observed cloud optical or 
microphysical properties (optical depth, albedo or cloud droplet effective radius) as a 
function of the change in the aerosol load (e.g. Twomey, 1974). When the number 

io concentration of aerosol particles capable of acting as cloud condensation nuclei in- 
creases, also the number concentration of cloud droplets is expected to increase. If 
the liquid water content of the cloud stays constant, the increase in the cloud droplet 
number concentration results in smaller average droplet size and thereby in a larger 
cloud optical depth and albedo. 

15 The aerosol-cloud interaction (ACI) is usually calculated using one of the following 
relations (Feingold et al., 2001): 


ACI= 
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d\na 


tflnr e 
LWP d\na 


LWP 


1 d\r\N d 
3 diner 


( 1 ) 


Here COT is the cloud optical depth, LWP is the cloud liquid water path, r e is the 
cloud droplet effective radius, N d is number concentration of cloud droplets, and a in 
20 some proxy for the aerosol burden such as the aerosol optical depth, total number 
concentration of aerosol particles or aerosol light scattering coefficient. The first two 
partial derivatives must be calculated at the constant LWP or cloud liquid water content. 

Aerosol number size distribution measurements give us an opportunity to investigate 
the effects of particle size on the aerosol burden and ACI. Therefore, when analysing 
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the ground-based data, or when combining ground-based in-situ measurements with 
retrieved cloud properties from remote sensing, we used the following equation for ACI: 


ACI (d k )=- 


d In r e 



_ 1 d lnA tot 

“3 / 30 \ ' 

d In I/V,- 

LWP \i=k / 


( 2 ) 


Here, A/, is the number concentration of particles in the DM PS size bin /, d k is the 
s particle diameter in the size bin k, the index k varies between 1 (c^sslOnm) and 28 
(c^s^S^nm), and /4 tot is the total number concentration of aerosol particles activated 
to cloud droplets. Equation (2) means that we calculate ACI by using the total aerosol 
number concentration in the diameter range d k to 500 nm as a proxy for the aerosol 
burden. 


io 3 Results and discussion 

Long-term average total particle number and sub-micron volume concentration at our 
measurement site are equal to 670cm -3 and 1.0 ^m -3 cm -3 , respectively (Dal Maso 
et al., 2008). This demonstrates that the site represents a very clean continental loca- 
tion. Measured cloud droplet number concentrations averaged between about 130 and 
is 1 80 cm -3 for the data sets 1 -4, which are at the lower end of values reported for conti- 
nental clouds (Miles et al., 2002). In the following, we will apply three different methods 
to estimate ACI at our site and then compare the results with available literature data. 

3.1 Ground-based data only 

Here, we used data sets 1-4 to investigate the behavior of ACI. To start with, let us 
20 consider a lower size limit of 50-1 50 nm when calculating aerosol-cloud interaction. In 
that case, ACI (d k ) varied in the range 0.2-0.3 depending on the exact value of d k and 
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used data set (Fig. la). The correlation between different data points used to calculate 
ACI was high (~0.9; Fig. 1 b), which is at least partly due to the fact that d k is close to 
the “dry” particle diameter that separates activated cloud droplets from cloud interstitial 
particles at our measurement site (Komppula et al., 2005). When these two diameters 
5 are close to each other, the particle population contributing to cloud droplet activation is 
roughly the same as the one used to calculate the aerosol burden. Values of ACI were 
clearly the lowest for the data set 3 in the range <4=50-1 5 nm. A major contributing 
factor for this feature was probably the different shape of the average particle number 
size distribution between the data set 3 and other data sets (Fig. 1 c). Data set three 3 
io had also the largest variation in the accumulation mode (Fig. Id), which may also 
contribute to higher ACI values at larger particle sizes. 

Above 150nm, ACI(<4) decreased rapidly with increasing d k , reaching values be- 
tween about 0.1 and 0.2 at d k of 370 nm (Fig. 1 a). The apparently lower values of ACI 
for <4>1 50 nm do not indicate weaker aerosol-cloud interaction, but rather demonstrate 
15 that the sub-population of aerosol particles used to calculate the aerosol burden is not 
anymore a good representative of the one that participates in cloud droplet activation. 
This view is supported by simultaneous strong decrease in the correlation between 
different data points used to calculate ACI (Fig. 1b, exception the data set 3). 

When <4 was decreased below 50 nm, values of ACI remained high, mostly 
20 above 0.25 (Fig. la). The main reason for this is probably that sub-50 nm particles 
did not give a big contribution to the total particle number concentration for the data 
sets 1-4 (Fig. 1c). Furthermore, the variability of the aerosol number size distribution 
was relatively low for these data sets (Fig. Id). 

3.2 Combined ground-based and satellite data 

25 For the data set 5 discussed here, values of ACI were calculated using either r e or COT 
obtained from satellite data as the cloud property. The aerosol burden was taken from 
ground-based particle number size distribution measurements in the same way as in 
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Sect. 3.1 . After filtering the data according to the criteria given in Sect. 2.1 .2, a total of 
62 cases were left for further analysis. 

The values of ACI for this data set were clearly lower than those obtained for data 
sets 1-4 (Fig. 2a), as were also the correlation coefficients between the different data 
5 points used to calculate ACI (Fig. 2b). This demonstrates that when aerosol and cloud 
properties are not exactly co-located, determination of aerosol-cloud interactions be- 
comes challenging. There was also a clear, almost systematic, difference of about 0.05 
in ACI when either r e or COT was used as cloud property. A difference between these 
two ways of calculating ACI has been observed also in a couple of earlier studies (Mc- 
io Comiskey et al., 2009; Pendithurai et al., 2009). 

Contrary to data sets 1-4, values of ACI in the data set 5 decreased below 0.15 
when the lower limit of particle sizes used to calculate the aerosol burden went below 
30 nm (Fig. 2a). The probable reason for this behaviour was the relatively large number 
concentration of particles smaller than 30 nm and especially its large variability (Fig. 2c 
is and d). By looking at individual particle number size distributions, particles in this 
size range could be traced to new-particle formation events, which is a common phe- 
nomenon around our measurement site during late summer and early fall (Dal Maso et 
al., 2007). This can be observed by high variation in size distributions (Fig. 2d). Prac- 
tically no new-particle formation events were observed in the data sets 1-4, since the 
20 weather conditions that cause the Sammaltunturi to be within a cloud are not generally 
favourable to atmospheric nucleation. 

3.3 Satellite data only 

Here, we relied solely on satellite data for determining ACI. The values of r e and COT 
were obtained as in Sect. 3.2, whereas the aerosol burden was replaced with the re- 
25 trieved aerosol optical depth (AOD). A total of 505 data points were obtained when 
using the constraints given in Sect. 2.1 .3. The data were then averaged into AOD bins 
of size 0.02(0-0.02, 0.02-0.04, 0.04-0.06 etc.). At least four data points in each bin 
were required to accept the bin to be taken into consideration. 
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Figure 3 shows the dependences of r e and COT on AOD. A clear aerosol indirect 
effect can be seen, with of r e decreasing and COT increasing for larger values of AOD. 
The value of ACI obtained from this data set was 0.1 for both r e and COT taken as a 
cloud property, the respective correlation coefficients being equal to 0.45 and 0.35. 

5 Interestingly, the value of ACI obtained for this data set is not very far from the values 
of ACI(cf fr ) in data sets 1-4 when taken 370 nm as the smallest size (d k ) to calculate 
the aerosol burden. The explanation for this can be found when looking at how opti- 
cal measures of the aerosol burden, such as the aerosol scattering coefficient, are re- 
lated to aerosol burden determined from particle number size distributions. Aaltonen et 
10 al. (2006) investigated three years of simultaneous aerosol scattering and particle num- 
ber size distribution data from our measurement site. They found that while the aerosol 
scattering coefficient had quite a good correlation with the total number concentration 
of accumulation mode particles (1 00-500 nm; R 2 =0.60), the respective correlation was 
much better with the number concentration of particles >200 nm in diameter (fl 2 =0.85) 
15 and even better with that of >300 nm particles (fl 2 =0.93). By putting these pieces 
of information together, we may conclude that, at least in our measurement site, using 
the aerosol scattering coefficient (and therefore also AOD) as a measure of the aerosol 
burden is very likely to lead to too low values of ACI. 

3.4 Comparison with other studies 

20 The values of aerosol-cloud interaction (ACI) for our in situ ground-based measure- 
ments were mostly above 0.25 when keeping the lower size limit to calculate the 
aerosol burden below 100-1 50 nm. McComiskey and Feingold (2008) summarized 
results from seven different in situ airborne studies made for determining aerosol-cloud 
interaction. In five of these studies observed values of ACI were comparable to ours 
25 (>0.20), whereas two other studies showed clearly lower values ( 0.1). 

By combining in situ ground-based and satellite data, and by leaving out the smallest 
particles to calculate the aerosol burden, we obtained ACI in the range of 0.10-0.14 
when using r e as a cloud property and in the range of 0.16-0.20 when using COT as a 
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cloud property. A few other studies have determined ACI with the help of ground-based 
remote sensing of clouds and aerosol scattering coefficient measurements. Those of 
them that used r e as a cloud property, reported values of ACI that were in the range of 
0.13-0.19 for Arctic clouds (Garret et al., 2004), in the ranges of 0.02-0.16 (Feingold 
5 et al., 2003) and 0.04-0.1 7 (Kim et al., 2008) for continental clouds, and in the range of 
0.02-0.18 for marine stratus clouds (Pandithurai et al., 2009). Corresponding ranges 
of ACI when using COT as a cloud property have been found as 0.05-0.16 for coastal 
marine clouds McComiskey et al. (2009) and 0.01-0.15 for marine stratus clouds (Pan- 
dithurai et al., 2009). Taken together, these measurements indicate that the combined 
io use of remote sensing and in situ data results in value of ACI that (i) have a relatively 
large variability even for the same environment, (ii) are almost always smaller than 
0.20, and (iii) are clearly lower than ACI obtained from in situ measurements. 

At our site, ACI determined using the satellite data was equal to 0.1 regardless of 
whether r e or COT was used as a cloud property. Other satellite-based investigations 
is have observed quite a variable range of values for ACI over different world regions (see 
Bulgin et al., 2008, and references therein). However, when averaging over the globe, 
Bulgin et al. (2008) found ACI to be in the range 0.10-0.16 during different seasons. 
This is line with our finding and suggests that satellite-derived values of ACI tend to be 
lower than those obtained using either solely by in situ data or a combination of in situ 
20 and remote sensing data. 


4 Conclusions 

In this work, we investigated the first indirect effect, or aerosol-cloud interaction (ACI), 
over a clean, northern high-latitude site using three different approaches. All of the 
approaches used were able to detect the first indirect effect at our site. The obtained 
25 values of ACI were highest for in situ ground-based data, clearly lower for combined 
ground-based and satellite data, and lowest for data relying solely on satellite retrievals. 
This feature is consistent with available literature data, even though it must be kept in 
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mind that observed values of ACI display usually a large variability even for a single 
site. 

Perhaps the most important of our findings was the high sensitivity of ACI to the 
definition of the aerosol burden. For our in situ data ACI was typically at its maxi- 
5 mum when the sub-population of aerosol particles used to calculate the aerosol bur- 
den was roughly the same as the one participating in cloud droplet activation, which 
meant counting all particles down to about lOOnm in diameter. Optical measures of 
the aerosol burden do not usually fulfil this requirement so, at least for our site, using 
the aerosol scattering coefficient or AOD as a measure of aerosol burden would result 
io in values of ACI that are smaller than the “real” aerosol-cloud interaction. 

Besides the aerosol burden, there are number of other factors that might affect the 
first indirect, or the variability of ACI as determined from measurements. One influenc- 
ing factor, observed also in this study, is that aerosol and cloud properties are often 
not measured at the same place or time (e.g. Myhre et al., 2007; Shao and Liu, 2009). 
15 It is also well known that the fraction of aerosol particles activating into cloud droplets 
depends on the air updraft velocity, shape of the particle size distribution and chemi- 
cal composition of the aerosol population (e.g. McFiggans et al., 2006; Reutter et al., 
2009). We observed a clear effect of the shape of the particle size distribution on the 
magnitude of ACI. 

20 McComiskey and Feingold (2008) made radiative forcing calculations and estimated 
that narrowing the uncertainty in measures of ACI to an accuracy of 0.05 would place 
estimated cloud radiative forcing on “a sounder footing”. This sounds challenging, es- 
pecially when given the large variability in reported values of ACI. We have shown that 
at least a part of the variability in ACI can be explained by how different investigators 
25 have related different cloud properties to “aerosol burden”. In this regard, more stud- 
ies should be dedicated to explore how well various optical measures of the aerosol 
burden can be used to calculate ACI in different environments. 
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0.35 




Fig. 1. (a) Aerosol-cloud interaction (ACI) in different cases as a function size where the sum of 
aerosol burden for ACI is taken from, (b) Correlation coefficient {R) as a function of size where 
the sum of aerosol burden for ACI is taken from, (c) Average particle number size distribution 
in different cases, (d) Standard deviation of the particle number size distribution in each case. 
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Fig. 2. (a) Aerosol-cloud interaction (ACI), as a function size where the sum of aerosol burden 
for ACI is taken from, (b) Correlation coefficient (fl) as a function of size where the sum of 
aerosol burden for ACI is taken from, (c) Average particle number size distribution, (d) Standard 
deviation of the particle number size distribution. 
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Fig. 3. On the left is the cloud is droplet effective radius (r e ) as a function of AOD. On the right 
is cloud optical thickness (COT) as a function of AOD. Data are based on satellite retrievals. 
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